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Abstract — Although  seismic  waves  have  been  studied  for 
many  years,  their  soliton  structure  has  only  recently  studied. 
Deformation  solitons  propagate  along  earthquake  faults  and 
can  induce  large  earthquakes.  Rotation  solitons  are 
generated  in  earthquake  sources  and  propagate  throughout 
the  Earth.  The  conclusion  to  be  reached  from  our  paper  is 
that  the  research  on  seismic  solitons  is  essential  for 
investigating  the  propagation  of  seismic  waves  and  helps 
understand  mechanisms  triggering  earthquakes.  This  paper 
discusses  the  development  of  elastodynamics  equations 
similar  to  Maxwell's  equations  in  a chiral  -mode  which  is 
applied  to  a seismic  channel,  which  is  dispersive  and 
nonlinear.  The  chirality  is  described  in  terms  of  the 
formalism  proposed  by  Born-Fedorov.  The  nonlinearity  is 
Kerr-type,  and  dispersion  of  the  medium  is  taken  into 
account  explicitly  through  the  Taylor  series  expansion. 
Through  numerical  calculations  these  theoretical  results 
allow  us  analyze  the  soliton  propagation  of  S-seismic  pulses 
which  can  induce  strong  earthquakes.  The  numerical 
calculation  is  applied  to  the  San  Ramon  Fault  localized  in 
Santiago  City,  Chile  which  is  a seismically  active  fault  that  is 
a main  element  to  be  considered  in  any  study  on  seismic 
hazard  assessment  for  this  city. 


Keywords — chirality,  S waves,  San  Ramon  Fault, 
earthquake,  soliton. 

I.  INTRODUCTION 

Paleoseismic  studies  of  a prominent  fault  scarp  at  the  west 
flank  of  the  Andes  in  Santiago,  Chile  (The  San  Ramon  Fault) 
were  published  recently.  The  San  Ramon  Fault  (SRF)  has 
been  recently  identified  like  a geologically  active  fault  which 
can  produce  large  magnitude  crustal  earthquakes,  in  the 
range  of  mw  U 6.9 -7.4  [1-6].  Figure  1 shows  the 

geomorphologic  map  of  the  San  Ramon  fault.  Seismic 
moment  can  be  estimated  with  a fixed  mpture  width  of  30 
km,  consistent  with  the  structure  of  the  SRF  at  depth  (Armijo 
et  al.,  2010;  Perez  et  ah,  2013)  and  5 m of  average  slip.  In 
that  case,  assuming  that  the  fault  scarp  is  either  15  km  long, 
according  to  the  most  emergent  part  of  the  scarp,  or  35  km 
long  if  considering  the  entire  well-recognized  SRF  (Armijo 
et  al..  2010).  predicts  Mw  U 7.25  and  Mw  0 7.5 , respectively. 

Those  magnitude  estimates  would  imply  surface  mpture  and 
significantly  stronger  ground  motion  (and  destmctive  effects) 
for  Santiago  (Perez  et  al.,  2013;  Pilz  et  al.,  2011). 


Using  5 m 

Fig.l:  The  San  Ramon  fault  scarp,  Chile  ( WAT — west  Andean  thrust).  The  geomorphologic  map  (Armijo  et  al.,  2010)  is 
draped  over  a 9 m digital  elevation  model.  The  San  Ramon  fault  forms  a continuous  scarp  that  limits  to  the  west  the  uplifted 
piedmont  covered  with  Quaternary  alluvium  and  incised  by  streams  (early-middle  Pleistocene  in  yellow  tones,  late 

Pleistocene-Holocene  in  white). 
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To  further  investigate  site  conditions  in  the  Santiago  basin 
we  must  consider  different  locations  of  possible  hypocentres 
along  the  San  Ramon  Fault,  a multikilometric  frontal  thrust 
at  the  western  front  of  the  Main  Cordillera,  which  has  been 
shown  to  pose  a significant  seismic  hazard  for  the  city 
(Armijo  et  al.  2010).  The  best  surface  expression  of  the  San 
Ramon  Fault  is  found  along  an  around  15  km  long  segment 
with  a sharp  fault  trace  in  the  25  km  long  part  separating  the 
rivers  Mapocho  and  Maipo  along  the  San  Ramon  mountain 
front  (see  Fig.  1). 

In  this  paper  we  propose  a model  to  estimate  the  magnitude 
of  earthquakes  along  the  fault  of  the  SRF  which  acts  as  a 
seismic  channel  which  develops  a soliton  pulse  acting  as 
earthquake  source.  Rotational  motions  in  earthquake  sources 
naturally  generate  rotational  seismic  waves  [7-10]. 

The  difference  in  speed  of  travel  of  P- waves  and  S-wave  is 
vital  to  transmit  energy  of  seismic  wave.  The  P wave  is  a 
longitudinal  wave  or  a compression  wave.  Force  is  applied  in 
the  direction  that  the  wave  is  travelling.  In  S wave,  the 
medium  is  displaced  in  a transverse  (up  and  down  - 
compared  to  the  line  of  travel)  way,  and  the  medium  must 
shear  or  "move  away"  from  the  material  right  next  to  it  to 
cause  the  shear  and  transmit  the  wave.  This  take  more  time, 
and  this  is  why  the  S wave  moves  more  slowly  than  the  P 
wave  in  seismic  events.  Here  we  model  the  SRF  like  a 
seismic  channel  using  a chiral  approach  [11-13]. 

There  are  two  types  of  tectonic  solitons  which  can  be  excited 
in  the  earthquake  source  and  propagating  along  the  fault: 
longitudinal  self-distortion  (plastic)  solitons  and  shear  self- 
distortion solitons.  The  Earth’s  interior  is  modeled  as  an 
elasto-plastic  continuum  and  elastoplastic  waves  [9,  10],  that 
describe  elastic  longitudinal  and  plastic  shear  distortion 
solitons.  We  apply  similar  method  in  order  to  determine  the 
solitons  propagating  along  the  fault  SRF.  The  type  of  soliton 
studied  here  is  seismic  shear  self-distortion  Schrodinger’s 
solitons  which  may  propagate  slowly  along  the  fault  and 
trigger  new  strong  earthquakes.  Our  main  results  (equation 
20)  is  obtained  under  the  chiral  approach  [11-13], 

The  chirality  was  first  observed  as  optical  activity  and 
corresponds  to  the  rotation  of  the  plane  of  polarization  in 
linear  isotropic  materials.  The  Born-Fedorov  equations  for 
electromagnetic  system  are  given  by  [14,  15]. 

D = e (e  + TcVxe)  (1) 

B = p(h  + TVxh)  (2) 

The  pseudoscalar  T represents  the  measure  of  chirality  and 

has  units  of  length.  The  validity  of  equations  (1)  and  (2)  has 
been  demonstrated  in  studies  of  chiral  electrodynamics  [11] 
and  from  an  electromagnetic  point  of  view,  chiral  soliton 


pulse  can  be  described  by  specific  equations  [12]-[13].  By 
analogy  between  elastodynamics  and  electrodynamics  in  this 
work  we  use  these  Born-Fedorov  equations  which  are  the 
most  suitable  to  study  SRF  like  a seismic  channel.  Here  the 
electric  permittivity  is  equivalent  density  of  matter £—>/?, 
the  magnetic  permeability  is  equivalent  to  the  inverse  of  the 
Lame  parameter  p — > 1/  jus . The  speed  of  transversal  seismic 

S-  waves  is  given  by  yjp/  ps  [14,  15]  and  changes  of  the 

polarization  of  transverse  seismic  waves  are  take  into 
account  during  their  propagation  through  the  SRF. 


II.  BASIC  EQUATION  FOR  ROTATIONAL 
PROPAGATION. 

Using  equations  (1)  and  (2),  we  obtain  in  the  analogous 
framework  under  this  section  the  nonlinear  Schrodinger 
equation  for  a chiral  seismic  channel.  (S-type  waves) 


Ds  = P„Es+pTc(VxEs) 

(3) 

Bs=l/ps(Hs+T\VxHs)) 

(4) 

These  equations  are  symmetric  under  duality  transformations 
and  temporal  reversibility,  pn  is  the  density  of  matter  and 


Tc  the  chiral  rotational  seismic  coefficient.  The 
corresponding  like  Maxwell  equations  are 


3 Ms) 


Vx//..  = 

5 dt 


<jE,+—T  (VxE,)=^^  + aEs+pTS/x^P- 

s 2,  CV  s>  2,  s r c 


dE, 


dt 

dH 


dt 


dt 


(5) 


dB  dH  d(Vx//  ) 

VxEs=--^-  = -(l/ps)-^-(l/ps)Tc  ^ (6) 

Taking  the  rotational  of  the  equation  (6)  and  considering  the 
following  approximations 

V-Ds  DO  = pty-Es+pTN  VxEs;  VTs=0->V/),  0 


V ■ £,,  = 0 =>  V ■ //  v = 0 
We  obtain  the  following  wave  equation 

V2Es  + (l/ps)pT2  = (1  lps)pn 

dE  d2E 

+(l/Ms)<r-^-+{(l/Ms)PnTc  + (l/Ms)P)TcVx-^f  (7) 

dE 

+(1/  ps  )oTV  x — A 
dt 

Here  we  assume  that  the  chiral  seismic  medium  is  of  a Kerr 
nonlinearity  type,  described  by  a refractive  index  such  that 
the  seismic  permittivity  is 

Pn  ~ Ps  Pi  \ES  | (8) 
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where  ps  is  the  linear  part  and  p , is  the  nonlinear 
component  respectively  of  pn  rx  is  the  rock  conductivity 
loss.  Past  studies  have  demonstrated  that  the  slip  type  or  style 
of  faulting  is  potentially  significant  for  correlating 
earthquake  magnitude  and  mpture  parameters.  To  categorize 
the  dominant  slip  type  for  SRF  in  this  paper,  we  use  a simple 
scheme  based  on  equation  (8). 

From  equation  (8)  we  can  be  inferred  easily  the  expression 
for  the  index  of  refraction  as  in  [8],  Replacing  pn  into 
equation  (7)  we  obtain 


-2  jk0—-k^ 
oz 


a1? 


2 jco0- — <an'P  = 


(11) 


32£ , 


d2E, 
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where  v2s  = ps  ^ ; a = (1  / ps  )cr ; kQv  = 0)Q ; 

P = (\l ps)p2 , j = yf— 1.  To  get  to  the  equation  (13)  the 

approximation  of  small  amplitudes  is  also  considered,  given 
by 
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Assuming  that  Es  represents  a localized  waveform,  which 
propagates  in  the  z direction,  it  has 

(10) 

where  lF  represents  the  complex  envelope.  To  solve  the 
equation  (9)  the  property  of  the  Fourier  transform 

3“ / dt2  <->  -a>~  is  applied,  and  then  we  determine  the 
operators  V"  and  V X . 

After  several  algebraic  manipulations  the  result  is  as  follows 


Es(r,t)  = (x+  jy)x¥(r,t)e~i{kz~%,)  = 


To  illustrate  the  difference  between  the  high  frequency  and 
the  low  frequency  in  an  earthquake  we  show  the  figure  2, 
which  corresponds  to  the  strong  earthquake  February  27, 
2010,  Chile,  ffere  we  have  a set  of  measurements  of  high- 
frequency  energy  radiation.  The  upper,  middle  and  lower 
traces  are  an  observed  seismogram,  the  squares  of  the  band- 
pass (2-4  Flz)  filtered  seismogram,  and  its  smoothed  time 
series  (normalized  by  the  maximum  value),  respectively.  “A” 
and  “F”  in  the  lower  trace  denote  the  arrival  of  P-wave  and 
estimated  end  of  high  frequency  energy  radiation, 
respectively  ( CO0  have  correspondence  with  the  upper  signal 

of  high  frequency>  10  radians/seg  and  — have 

3 1 

correspondence  with  the  lower  signal  of  slow  frequency)  [14, 
18].  In  our  analysis  we  assume  that  a possible  earthquake  in 
San  Ramon  Fault  waveforms  can  be  similar  but  of  lower 
amplitude  to  those  shown  in  Figure  2 


Fig. 2:  A measurement  of  high  frequency  energy  radiation  of  the  February  27,  2010  Chile  Earthquake,  the  bottom  graph  shows 
the  low  frequency  signal  which  is  the  envelope  of  the  high  frequency  wave.  In  the  next  section  this  envelope  is  modeled  as  a 

Gaussian  pulse. 
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The  effect  of  dispersion  is  included  heuristically.  Making  the 


change  of  variables  (j)  = 2k  y Z =■ 


and 


rearranging  the  terms  of  equation  (11)  we  obtain 

j 


d</>  1 d<j) 

dz*  v dt 


= -(l -k0Tc)^<f> 


2k 


M2 

(2  k): 


(i-kjMt 


(12) 


As  the  envelope  VT  ( j;,  i ) is  a slowly  varying  function  of  in 
Z and  t,  the  dispersion  relation  k=k[oi)  can  be 
transformed  to  the  domain  of  spatial  variations  by  means  of 
A (0=  CO—  (O0,  which  is  a small  deviation  of  the  sideband 

frequency  with  respect  to  0)Q  , and  through  A k = k —kQ, 
which  represent  the  corresponding  wave  number. 

Using  the  Fourier  transform  for  A k , A (0  , approximating 
A 2 -ji-  5 and  using  the  Taylor  series  we  obtain 


13  dk  3 . 1 3 A:  3 

v 3 1 dot  dt  ^ 2 3 of  dt 2 


ai  dt 


By  substituting  this  operator  in  equation  ( 14),  we  obtain 
dz  vg  dt  2 dt  ' ° c>  2 k/  K ° c\2 kf 

+(i-^Wo 


(13) 


where 

k'-  — - — - k"~ 

3 O)  vg  ’ dof  ’ 

Equation  (13)  describes  the  propagation  of  pulses  in  a chiral, 
rotational,  dispersive  and  non  linear  channel  [13]. 

The  first  term  represents  the  evolution  of  the  pulse  with 
distance.  The  second,  and  third  terms  represents  the 

scattering  of  a seismic  chiral  channel  &r|=-^-jand  k" 

corresponds  to  the  chromatic  dispersion  k' , indicating  that 
the  pulse  moves  with  the  group  velocity,  while  the  dispersion 
of  the  group  velocity  (GVD)  is  represented  by  k " , which 
alters  the  relative  phases  of  the  frequency  components  of 
pulse  widening,  producing  its  temporal  expansion.  The 
fourth  is  associated  with  attenuation  ( OC ),  in  this  case  such 
losses  are  weighted  by  the  seismic  chirality.  The  fifth  term 

\ff)\  0 represents  nonlinear  effect,  like  a Kerr  effect,  which  is 


characterized  by  having  a refractive  index  dependent  on  the 
seismic  field  intensity  and  related  to  the  slip  type  that  can 
alter  and  widen  the  spectrum  of  frequency.  This  term  also 
depends  on  the  chiral  factor.  Finally,  the  last  term  is  clearly 
associated  with  the  chirality  of  the  seismic  channel. 

III.  ANALYSIS  AND  SIMULATIONS 

In  order  to  ease  up  the  solution  of  the  propagation  equation 

Z 

the  following  change  of  variables  is  introduced:  t*  = t , 

vg 

z*  = Z 

Defining 

v2s=Ms± ; a = (\//is)a-.  k0  = ^-,  J3  = (1/ jus)p2 


q = 2£q(Z>  , Z = 


C = l~k0Tc 


(14) 


oi. 
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i -CK 
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k - ■ 
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d2k 
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k dk  / d(0 


(15) 


We  have 
3 q 


J — = — k 
3 Z 2 


a2 
d q 

dt*2 


co^aC 

2 k r, 


q + 


fko^a 
(2  kQ)3 
(16) 


<7  <7 


Finally,  is  useful  normalize  the  Eq.  (16)  by  introducing 


t* 

t = — 
tn 


(17) 


Let  us  model  heuristically  the  relationship  between  the  wave 

amplitude  and  the  initial  seismic  power  as  follows 

The  normalized  amplitude  q (equation  16)  is  proportional  to 


P . \P  4 f 

ref  ) ref  ref 


1 


Jref 


Where  Ps  is  the  power  earthquake,  Pref  is  the  power  reference 


and  t is  the  time  duration  of  the  earthquake,  / 


lref 


using  the  empirical  equation  of  Gutemberg  and  Richer  [16] 
which  is  useful  for  magnitude  of  shallow,  small-moderate 
nearby  earthquakes 
log  Es  = 11.8  + 1.5M 

(18) 

We  have 

Es  _ 

Eref 
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We  use  this  empirical  equation  because  this  theory  gives  a 
relationship  between  the  magnitude  Ms  and  the  energy 


released  Es  which  only  tends  to  saturate  for  really  great 
earthquakes. 

Thus  we  can  model  the  S-wave  q through  the  transformation 

/o.12(Ms-M  ) 2 kn 

q(Z,T)  = e'1  4 e 0 U(Z,z) 

(19) 


where  Ms  is  the  peak  magnitude  of  the  soliton  and  M is 

the  reference  magnitude  for  a weak  seism  which  acts  as 
calibration  standard.  Finally  the  expression  (16)  can  be 
scaled  as 

2 2 i .-v 

dll  1 k"  d~U  /3d)  lo.n (Ms-Mr.)  kn  \ |2 

j — = + he  e \u\  u 

dZ  2 t0  dz  (2 k0) 

(20) 

To  our  knowledge  this  is  the  first  time  that  the  magnitude 
M s of  earthquake  is  considered  into  the  nonlinear 

Schrodinger  equation.  Through  numerical  calculations  these 
theoretical  results  (eqs  14-20),  would  allow  analyze  the 
effects  of  chirality  on  the  attenuation  of  equation  (20)  for 
propagation  of  S-seismic  pulses.  Also  we  take  into  account 

the  magnitude  Ms  of  moderate  earthquakes  along  the  SRF 

which  can  trigger  strong  earthquakes  in  Santiago  City. 

To  solve  the  pulse  propagation  in  nonlinear  dispersive  media 
we  use  the  split-step  Fourier  method.  The  relative  speed  of 
this  method  compared  with  most  FD  methods  can  be 
attributed  in  part  to  the  use  of  finite-Fourier-transform  (FFT) 
algorithm  [17].  In  general,  linear  and  nonlinear  parts  of 
equation  (19)  act  together  along  the  length  of  the  earthquake. 
The  split-step  method  obtains  an  approximate  solution  by 
assuming  that  in  propagation  the  seismic  field  over  a small 
distance  Ac , is  carried  out  in  two  steps.  In  the  first  step,  from 
Z = Z0  to  Z = Z0  + AZ/2 , the  linear  part  acts  alone,  and  the 
non  linear  part  is  zero.  In  the  second  step,  the  nonlinearity 
acts  alone  in  the  point  AZ/2  , (linear  part  is  null). 

For  the  simulation  we  consider  an  average  of  velocity  vs 
because  the  bedrock  of  SRF  is  layered  with  vs  = 2.500ml  s 
between  0 and  2.2  km  depth,  vs  = 3.300m/ s between  2.2 
km  and  8 km  depth  (Barrientos  et  al.  2004)  [19].  The 
frequency  w0  was  set  to  > 10 radians / s in  this  simulation. 

The  wave  comes  to  8 km  deep  in  the  bedrock. 

Figure  3.  shows  the  Intensity  U calculated  numerically  from 
equation  (20).  The  parameters  correspond  to  a weak 
earthquake  with  m s = 1 , Mref  is  the  reference  magnitude  for  a 


weak  seism  which  acts  as  calibration  standard  to  nonlinear 

Schrodinger  equation  with  Ms  > 1 .Flere  rc  = 0 , C = 1 the 

chiral  effect  is  not  considered.  The  angle  0 between  the 

plane  which  contains  U and  axis  Z,  0 = 0 .The  normalized 

attenuation  is  neglected  for  short  distance,  Z □ 30 km . 
K0 

The  equation  (20)  is  reduced  to 

dU  1 k"  d2U  /3a>:  | ,2 

u 

dZ  2 tg  dr  (fkg) 

Which  is  the  typical  result  obtained  by  other  authors  [7,  9]. 


Fig.  3:  Parameters  of  a weak  earthquake;  M =1,  M s = 1 , 
Vs  =W0’Tc  =0’  C = l 

yjpl  jus  = 3 km  I s ps=4g/cmA3  p,  = 10~4  Z0  = 30km , 

J3  = (1/  ps)p2  = lcf 3 , C = 1 , 0 = 0,  depth  = 8km 

The  propagating  shear  wave  is  represented  by  a Gaussian 
pulse  simulating  the  bottom  graph  of  figure  2 
Figure.4.  shows  the  envelope  of  U for  a strong  earthquake 
with  C = 1 - k0Tc  = 0.3 , nearly  of  Z = 30 km  we  have  a peak  of 

U,  which  it  indicates  that  the  maximum  of  U can  trigger  a 
major  earthquake.  From  equation  (20)  we  can  obtain 

M s ~ M rtt  ~ 6-4 . that  is  M s ~ 7.4  . This  value  agrees  with  the 

estimations  of  Vargas,  Armijo  y Perez  [1,  2,  3]  respectively 
who  are  used  the  Wells  - Coppersmith’s  semi  empirical 
approach  [6].  We  note  that  in  this  paper  mw  □ Ms  because 

M s only  tends  to  saturate  for  really  great  earthquakes 
Mw  > 8.0  . The  properties  of  rock  in  the  seismic  wave  travel 

path  are  important  because  they  directly  influence  the 
amplitudes  of  seismic  waves.  Rock  is  normally  considered  to 
be  stable  with  respect  to  seismic  amplification,  but  here  we 
show  that  this  is  not  always  the  case.  The  wavelength  of  a 
seismic  wave  of  a given  frequency  is  proportional  to  its 
velocity,  v = vs  = Af  . If  a shear  wave  passes  into  a lower 
velocity  medium,  it  will  reduce  its  wavelength  so  to  conserve 
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energy,  the  wave  is  amplified.  Sedimentary  basins  have 
lower  velocities  than  bedrock  and  can  therefore  produce 
intense  amplification  of  incoming  seismic  waves.  As  seismic 
wave  travel  faster  through  harder  denser  rocks  and  compact 
sediments  than  through  soft  rocks  and  loose  sediments. The 
change  in  velocity  also  changes  the  rate  at  which  seismic 
energy  is  amplified.  As  a wave  travels  up  to  the  layer  with 
sediment  (nearly  of  Z = 30 km  ),  density  decreases  (as  does 
the  impedance),  therefore  the  shear  wave  velocity  also 
decreases.  If  Energy  is  to  remain  constant  when  wave 
approaches  the  surface,  the  particles  velocity  must  increase 
to  compensate  the  decrease  in  density  and  shear  wave 
velocity  producing  a large  effect  on  the  basin  layer.  Sediment 
impedance  affects  the  amplitude  of  earthquake  ground 
motion.  When  Z > 30 km , the  S-wave  along  the  rock  layer 
decreases  due  to  the  drastic  change  of  impedance  between 
the  rock  and  the  sediment  layer.  Impedance  is  the  product  of 
the  density  ( ps ),  the  shear  wave  velocity  ( vs ) and  the 
cosine  of  the  angle  of  incidence  which  is  defined  as  the  angle 
between  the  vertical  and  the  direction  of  seismic  wave 
propagation.  In  terms  of  angle  6S  the  change  in  impedance  is 
expressed  with  the  impedance  contrast  Ic  = psVs  sin  ds  / phVh 
At  Z0  = 30km  the  reflection  coefficient  is 
rc  = Psvs  sin  es  - PX  /(PsVs  sin  6S  + pbVb  ) with 
Ic=psVssm0s/pbVb»l  the  S-wave  has  a peak  with 
M s > 7 , that  is  the  energy  E reaches  a maximum  that  is 

transferred  to  the  basin  layer  where  seismic  waves  are 
trapped  and  begin  to  reverberate. 

Thus  we  can  see  that  the  chiral  factor  C contributes  to  the  S 
wave  to  be  amplified.  Consequently  the  S wave  travels 
further  along  the  fault. 


Fig.  4:  Parameters  of  a strong  earthquake,  Peak  of  M s > 7 , 
M,r,  ='■  C = l~k0Tc  =0.3 

<Jpl ps  = 3km/ s , ps  = 4g /cm  A 3 , p,  = 10(— 4)  Z0  = 30 km , 
w > 10  radians/s.  0 = 15° , depth  = 8knr 
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As  we  can  see,  rotational  tectonic  waves  can  propagate  along 
the  fractured  tectonic  fault  like  SFR  with  a speed  of  2-4  km 
per  second.  These  waves  may  have  a form  of  rotational 
seismic  solitons  and  can  trigger  major  earthquakes.  Thus,  the 
research  on  rotational  seismic  solitons  is  essential  for 
investigating  the  propagation  of  seismic  waves  and  helps 
understand  mechanisms  triggering  earthquakes.  Chiral 
rotational  seismic  waves  propagate  faster  in  solid  rocks  and 
much  slower  in  fractured  media  along  tectonic  faults.  A 
rupture  breaking  the  free  surface  is  another  important 
element  to  keep  in  mind  in  order  to  improve  future  kinematic 
rupture  scenarios  for  the  SRF  region.  Indeed.  In  this  model 
further  improvements  should  include  local  site/soil  effects 
and  basin/topographic  effects. 

IV.  CONCLUSIONS 

In  this  paper  we  have  obtained  the  Schrodinger  nonlinear 
equation  for  a channel  whose  core  is  chiral,  dispersive  and 
has  nonlinear  behavior.  The  effect  of  chirality  is  manifested 
by  the  terms  associated  with  the  loss  of  the  chiral  channel 
and  the  nonlinear  coefficient.  The  most  significant  result  of 
our  work  is  that  to  use  the  chirality  of  S waves  can  cancel  out 
the  losses  and  nonlinearities  of  the  channel,  which  would 
allow  us  to  modify  radically  their  behavior  as  channel  of 
seismic  soliton.  Numerical  calculation  of  the  equation  was 
obtained,  to  model  by  soliton  in  the  SRF.  We  propose  that 
the  San  Ramon  Fault  is  not  only  a geologically  active  fault, 
but  also  it  presents  seismicity  that  can  be  associated  with  this 
structure  at  and  resulting  in  a seismically  active  fault  too.  As 
seismic  waves  travel  though  the  ground,  they  travel  faster 
through  hard  rock  than  soft  soil.  As  a result,  when  the  waves 
move  from  hard  rock  to  soft  soil,  the  amplitude  (largeness)  of 
the  waves  needs  to  increase  to  be  able  to  carry  the  same 
amount  of  energy,  creating  stronger  shaking.  This  same 
principle  accounts  for  the  site  effects  of  sediment  thickness. 
The  deeper  the  sediment  above  bedrock,  the  more  soft  soil 
there  is  for  seismic  waves  to  travel  through,  therefore 
creating  stronger  amplifications.  This  is  a main  element  to  be 
considered  in  any  study  on  seismic  hazard  assessment  for  the 
city  of  Santiago,  Chile.  Through  numerical  calculations  our 
theoretical  results  would  allow  analyze  the  propagation  of  S- 
seismic  pulses  that  can  trigger  strong  earthquakes  in  this  city. 
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